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Abstract. Single-particle dynamics of the Anderson impurity model in the presence 
of a magnetic field H are considered, using a recently developed local moment approach 
that encompasses all energy scales, field and interaction strengths. For strong coupling 
in particular, the Kondo scaling regime is recovered. Here the frequency (w/wr) 
and field (H/lok) dependence of the resultant universal scaling spectrum is obtained 
in large part analytically, and the field-induced destruction of the Kondo resonance 
investigated. The scaling spectrum is found to exhibit the slow logarithmic tails 
recently shown to dominate the zero-field scaling spectrum. At the opposite extreme 
of the Fermi level, it gives asymptotically exact agreement with results for statics 
known from the Bethe ansatz. Good agreement is also found with the frequency and 
field-dependence of recent numerical renormalization group calculations. Differential 
conductance experiments on quantum dots in the presence of a magnetic field are 
likewise considered; and appear to be well accounted for by the theory. Some new 
exact results for the problem are also established. 
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1. Introduction 

The Anderson impurity model (AIM) [|IJ has long occupied a central role in condensed 
matter theory. Reviewed comprehensively in ||, it serves as a paradigm for the 
physics of strong local interactions, and remains the canonical model for understanding 
magnetic impurities in metals; competition between on-site Coulomb repulsion and 
band hybridization generating the Kondo effect in strong coupling, where the former 
dominates the latter. Renewed interest in the problem has arisen recently from the 
discovery that direct mesoscopic realizations of AIMs may be tailor made: quantum 
dots H ^| for example, or surface atoms probed by scanning tunneling microscopy 
Moreover the intrinsic tunability of such nanoscale devices offers controlled access to a 
wider range of quantum 'impurity' physics than usually accessible with more traditional 
materials. 
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Experimental probes of quantum dots are of course typically dynamical; and therein 
lies a major theoretical challenge. Static (thermodynamic and related) properties of 
the AIM are well understood, using a variety of powerful techniques that include the 
numerical renormalization group (NRG) ||, Fermi liquid theory 0, and the Bethe 
ansatz ||. That is not however the case for dynamical properties, such as the single- 
particle excitations. The many theoretical approaches to such are approximate by 
necessity, and suffer from well known qualitative limitations f|| - even for the AIM 
in equilibrium, let alone non-equilibrium effects that, to a greater or lesser extent, 
are relevant to e.g. non-linear differential conductance measurements || |J. And 
the fact that numerical methods, notably the NRG, can now provide benchmark 
numerical results for equilibrium dynamics (see e.g. @) gives added impetus to further 
development of approximate many-body theories. 

In this paper we consider the AIM in the presence of an applied magnetic field, H. 
This is a topical issue [9-15], motivated in part by differential conductance measurements 
on quantum dots in the Kondo regime |3], |4j] : these show a characteristic field-induced 
splitting of the Kondo resonance ||, indicative of the slow crossover from a locked Kondo 
singlet to an asymptotically free local moment with increasing h = We pursue 

it here within the framework of the recently developed local moment approach (LMA) 
[16-19], a technically straightforward non-perturbative many-body method in which the 
notion of local moments |I[] is introduced explicitly and self-consistently from the outset. 

The LMA handles single-particle dynamics on all energy scales, and for all 
interaction strengths U = U / nA |16], [HJ (with U the on-site Coulomb repulsion and A 
the hybridization strength). Most importantly, in strong coupling U 3> 1 it captures the 
Kondo or spin fluctuation regime characterized by a low-energy Kondo scale c^k that 
is exponentially small, such that zero-field single particle dynamics exhibit universal 
scaling in terms of uj/uj-^. The resultant scaling spectrum for the symmetric AIM yields 
good quantitative agreement with NRG results |2(], recovering not only Fermi liquid 
behaviour on the lowest energy scales, but also revealing slow logarithmic tails that in 
fact dominate the scaling spectrum [18]. 

For h ^ 0, LMA results for static properties have recently been considered 
|l"9|l , notably for the impurity magnetization and corresponding /i-dependent spin 



susceptibility in the Kondo regime. These too are found to yield good agreement with 
exact results known from the Bethe ansatz || |2~lH , being asymptotically exact in both 
the weak- and strong-field limits; and correctly recovering the field-independence of the 



Wilson ratio, R^(h) = 2 for all h ]S|, 32|. Here by contrast we consider field-dependent 
spectral dynamics, beginning (§2) with a brief introduction to the LMA for finite h; 
and focussing on the symmetric model, in which the spectral effects of an applied field 
are most simply apparent. There are two essential domains of field strength. First the 
important strong coupling, universal Kondo limit, where h by definition is irrelevantly 
small compared to the 'bare' electronic scales A or U; but with h/u^ arbitrary, and 
spanning the entire range of field strengths relevant to the Kondo model. It is this regime 
that is of primary interest in relation to experiments on quantum dots ||, Q . There are 
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also non-universal regimes of field strength, including a striking spectral signature of the 
crossover to the high-field limit where the fermions become effectively spinless. These 
are discussed as part of §3, where evolution of the AIM spectra on all frequency and field 
scales is considered, with the aim of extracting a broad overall picture of the problem. 

In §'s 4,5 we turn to the strong coupling Kondo limit, and spectral scaling in terms 
of uj/uk and h/u^. As for the zero- field case considered recently [I3|, our initial aim 
(§4) is to obtain analytically the /i-dependent scaling spectrum in a manner that is 
largely independent of the details of the LMA; and in particular to deduce explicitly 
(§4.1) the behaviour of the high-frequency spectral tails, as well as the field- dependence 
of the spectrum at the Fermi level. The main body of results for the LMA scaling 
spectra are given in §5; including comparison to recent NRG |14] and density matrix 
NRG (DM-NRG) (15], ^| calculations, as well as to the spinon approximation in which 



11| the single-particle spectrum is approximated by the density of states for spinon 



excitations obtained from the Bethe ansatz. The LMA results of §5.2 are found to 
compare favourably with the former, but are markedly at odds with those of the spinon 
approximation. The latter is discussed explicitly in §5.3, and its qualitative limitations 
identified. 

Finally, comparison is made in §5.4 to recent differential conductance experiments 
on quantum dots in a magnetic field; which, notwithstanding the natural absence 
of non-equilibrium effects in the LMA, are found to be rather well described by the 
theory. A brief summary is given (§6), together with an appendix in which, using 
microscopic Fermi liquid theory, we obtain two exact results for the problem that to our 
knowledge are new; specifically for the field-dependence of the quasiparticle weight, and 
the asymptotic low-field behaviour of the spectral shifts. 



2. Background 

The Hamiltonian for the AIM Q is given in conventional notation by 

H = ^ + ^ (eia + ^fli-a J fl i(J + ^ V ik (c\ a C ka + ll.C. j . (2.1) 

The first term refers to the host band of non-interacting electrons, with dispersion 
efc. The second refers to the impurity, with on-site interaction U and site-energy 
tin = £i — oh. The latter includes a local Zeeman coupling to the external field H 
(applied for convenience in the — z-direction), with h = ^gfi^H and o = +/— for 
t / |-spin electrons. The final term in equation (|2.1|) is the one-electron host-impurity 
coupling. For the symmetric AIM that we consider, = — ~, and by particle-hole (p-h) 
symmetry rij = ^2 a {ni a ) = 1 for all U and h. 

We focus on the total impurity Green function G(u; h) (with corresponding spectral 
density D(u\ h) = — 7r _1 sgn(a;) lmG(u; h))] where 
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with G a (u; h) = G^(u; h) — i sgn(u;)7rD (T (u;; h), given by 



G a (uj, h) = u + — A{uj) + ah — E CT (a;; h) 



-i 



(2.3) 



and uj + = L) + iO + sgn(co>). Here A(u;) = Ar(w) — i sgn(co>)Ai(c<j) (= — A(— u)) is the 
host-impurity hybridization, with Ai(a>) = ^Ylk^ik^^ ~ e k)i an d the hybridization 
strength A = Ai(u; = 0) is thus defined, with uj = the Fermi level. tl a (u;h) = 
S^(c<j; h) — i sgn(u;)Sj r (u;; h) denotes the impurity self-energy (excluding the trivial 
Hartree contribution that precisely cancels £j = — ¥)■ By p-h symmetry 



S CT (^; h) = -XL^-a;; h) 



[2A) 



and likewise D a (uj;h) = D_ a .(—u;h); whence D(u;h) = ^^2 a D a (uj;h) = D(—u;h) is 
naturally symmetric in u about the Fermi level. The /i-dependent quasiparticle weight 
Z(h) is defined by 

-l 



Z(h)= l-{d^{u J -h)/duj) ul=Q 



;2.s) 



(and from equation (|2.4j ) is a- independent). 

In practice we consider explicitly the usual wide-band AIM for which Ai(u) = 
A V uj, and Ar(u;) = 0. This is not of course restrictive since in strong coupling, 
U = U/tcAq ^> 1, the relevant low-energy Kondo resonance is a universal function 
of lj/uk (with the Kondo scale u>k defined as the HWHM of the h — Kondo 
resonance). It is thus independent of the detailed one-electron structure of A(o>)/A , 
which affects only the dependence of uk itself on the bare one-electron parameters 
(such as the bandwidth, D, of Ai(u)). By the same token the /i-dependence of the 
Kondo resonance is independent of whether a local or global magnetic field is considered. 
Application of a global uniform field leads additionally to an /i-dependent hybridization 
A a (u; h) = A(u> + ah); but its /i-dependence arises on a one-electron energy scale, and 
is thus irrelevant in the Kondo regime of finite h/uja and wk oc exp(— 7r?7/8Ao) — > 0. 

Within the LMA [16-19] the self-energy T, a (uj; h) is separated as 



t a (uj; h) 



a 



U\v(h)\ + X a (uj;h) 



(2.6) 



into a static Fock contribution (with local moment \fi(h)\) that alone would survive at the 
simple mean-field (MF) level of unrestricted Hartree-Fock; together with a dynamical 
contribution H a (w,h). The latter includes in particular a non-perturbative class of 
diagrams (figure [I]) that embody dynamical coupling of single-particle excitations to 
low-energy transverse spin fluctuations. Other classes of diagrams may also be included, 
but retention of the dynamical spin-flip scattering processes is essential to capture the 
strong coupling Kondo regime for U ^> 1 [|T^, [H]] . These are expressed in terms of MF 
propagators (solid lines in figure |l|) , viz 



Gcr(u\ h) 



u + -A(uj)+a (-U\n(h)\ + h 



-i 



[2.7) 
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Figure 1. Principal contribution to the LMA E cr (a;), see text. Wavy lines denote U. 



with spectral densities D®(u) = —n 1 sgn(u;) lmQ a (u; h); and £f(o;; fo) 
is given explicitly by [ p~6f| 

oo 



-uj; h)) 



; *0] (2. 



with Q^{u\h) the one-sided Hilbert transforms of Q a {oj\K) and 0(a;) the unit step 
function. Here, IT (&>;/&) denotes the transverse spin polarization propagator (shown 
hatched in figure ||). It is given at the simplest level by an RPA-like particle-hole ladder 
sum in the transverse spin channel, i.e. IT = °IT / (1 — U°U + ~) with °n + ~(c<;; h) the 
bare p-h bubble, itself expressed in terms of MF propagators. 

The central physical idea behind the LMA for h = is that of symmetry restoration, 
as detailed in [IB, IT71 : restoration of the broken symmetry endemic at pure MF level, 
via the spin-flip dynamics embodied in E cr (c<j;0). The MF propagators form the basis 
for constructing the dynamical self-energies = So-^] shown in figure [I]; and by 
themselves correspond to local symmetry breaking for |//(0)| > 0. Symmetry restoration 
is embodied mathematically in E?(0; 0) = £|"(0; 0), i.e. using p-h symmetry by 

1 



(0;0) 



U\n(0)\ + Sf(0;0) =0. 



(2.9) 



Equation ( |2.9|) is imposed self-consistently, achieved in practice |L6|, |T^| for given 
U by varying the local moment |/x(0)| from its pure mean-field value | /xo | (where 
yUo — (wjf — fiii)o with the average over the MF ground state). Symmetry restoration 
ensures correctly that the \ow-u behaviour of G(u; 0) constitutes a renormalization of 



the non-interacting limit; as manifest [I6| in a resultant quasiparticle form for G(uj; 0), 
the recovery of Fermi liquid behaviour at low energies, and preservation of the U- 
independent pinning of the Fermi level spectrum (viz ttA D(0;0) — 1 V U > 0, as 
follows directly from equations (P72L |2.3p using equation (|2.9|)). 
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Most importantly, self-consistent imposition of equation ( |2.9| ) introduces naturally 
a low-energy spin-flip scale u;^ = u) m (h = 0), evident in particular in a strong resonance 
centred on w = in the transverse spin polarization propagator Imll + ~(ct;; 0). Its 
form in strong coupling is readily deduced from equation ( |2.8| ) for u = = h, which for 
U ^> 1 has the asymptotic form [[lj] 



Sf(0;0) = ^ln 

7T 



A 



(2-10) 



with A = min[~,Z)]. Combining equation (j2.10 ) with the symmetry restoration 



condition equation (|2.9| ) (using |/i(0)| — > 1 in strong coupling) gives ~ 
A exp(— 7tU/8Aq). This is the Kondo scale, exponentially small in strong coupling and 
recovering the exact exponent 0. It naturally has no counterpart at simple MF level, 
and within the LMA the physical significance of the Kondo scale u;^ (oc co>k) is that it 
sets the timescale (~ fi/wj,) for restoration of the broken symmetry inherent at pure 
MF level. 

Before turning to h ^ 0, note that for h = two degenerate MF states arise 
TEL reflecting simply the invariance of H under a — > —a (f / |-spin symmetry). 



These states, denoted temporarily by a = A or B, correspond respectively to a local 
moment of fi = ±|/x|. The corresponding MF propagators are then denoted in full 
by Gaa{u;0), the self-energies constructed from them by Ti aa (u>] 0), and the resultant 
many-body Green functions by G acr (uj;0); for example, equation ( |2.2j ) in full notation 
is G(u>; 0) = \ Yl G aa (u>; 0) with the G afT (uj;0) given from equation (|2.3| ) in terms of 

Il aa (uj; 0). In writing equations ( |2.2[ - |2.5| ) we have assumed implicitly that either of the 
two mean-field states may be employed for h = 0. This is indeed correct, for the 
G aa {uj] 0)'s are related by ] / |-spin symmetry, viz 



Hence 



G Act (lj;0) = G B -*(u;0). (2.11) 



G(w;0) = i^G? a<T (a;;0) (2.12) 



is independent of a (which is not moreover specific to the symmetric AIM, since p-h 
symmetry has not been used in any way). Using equation ( [2.11| ), G(u; 0) may also be 
written equivalently as 

G(c;0) = ^G QCT (^0). (2.13) 

a 

This form shows that G(u; 0) may be viewed equivalently as involving an average over 
the two degenerate MF states a = A, B. Equation ( |2.13| ) is moreover independent 
of spin, a, reflecting the fact that G(uj\ 0) is equivalently the cr-spin Green function. 
Hence G(u; 0) = [u + — — £(c<j; 0)] _1 in terms of the conventional single self-energy 
H(u;0) which, using equation (|2.13|) , may therefore be obtained from the underlying 



two-self-energy description inherent to the LMA as discussed in [16-18]. 
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The above situation naturally changes for h 7^ 0. The degeneracy is removed and 
one or other MF state is picked out, according to sgn(ft). G(u; h) is then given by 

G(w,h) = ^Y,G cu ,(u;h) (2.14) 

a 

with a = A for h > and B for h < (given our convention for the Zeeman coupling 
in the Hamiltonian, equation |2J]). But the invariance of H under (cr, h) — > (—a,—h) 
implies 

G ACT (^|/i|) = GB- ff (^-|/i|) (2.15) 

(and likewise for the T, a(T {uj;h)), and thus G(u>;h) = G(u>;—h). Only ft, > need 
therefore be considered, and hence a = A; this will be assumed henceforth (and is 
already implicit in equations (|2.2| - [2.10|) ), and the a label thus dropped. Finally, note 
that while G(u;h) evolves continuously in h to its h = limit G(u;0), the fact that 
the degeneracy of the MF states is removed for h ^ means that G(u; h) for h ^ is 
naturally not expressible in the form of equation ( |2.13| ). 



3. Dynamics: all scales. 



We first consider single-particle dynamics on all frequency and field scales, i.e. 
encompassing both the low-energy Kondo resonance and high-energy Hubbard satellites, 
as well as the full range of magnetic field strengths. The /i-dependence of the low-energy 
scaling spectrum that arises in the Kondo limit will be pursued in the following sections. 

G(u>; h) is given by equations (|2.2|) , and using equations ( |2~3| , |2l| [2~7D G a (uj; h) may 
be expressed as G a (u>; h) = [Q~ x (uj; h) — T, a (uj; h)]^ 1 in terms off the dynamical part of 
the self-energy, S CT (o;; h). The latter, given in practice by equation ( pT8|) and figure [l], is 
quite generally a functional of the MF propagators Q a {oj\li) given by equation (|2.7|) : 
Ti a {uj\ h) = T: a [Q a ]. Hence G a (uj;h) = G a [Q a ], and from equation ( |2.7| ) it follows 
that the ft-dependence is embodied fully in x{h) = ^U\fi(h) \ + h. The LMA G{u;h) 
for h > is thus formally equivalent to that for h = 0, but with x(0) replaced by 
x{h) = x(0) + [^U5\fi(h)\ + h}; where S\fi(h)\ = \^(h)\ — |/i(0)|, and the zero-field moment 
|/i(0)| is determined such that symmetry restoration (equation (|279|) ) is satisfied. 

In practice we take 5\fi(h)\ to be given approximately at MF level, specifically for 
the wide-band AIM by 



7T 



5\pi(h)\ ~ - tan" 1 - tan- . (3.1) 



A J \ A J 



The ft-dependence of 5\fi(h)\ is in fact of little importance, since it is never significant 
compared to the bare Zeeman term (h) in x{h) — x(0) = ^U6\fi(h)\ + h. In strong 
coupling (SC) U = U/nA ^> 1, this is seen directly from equation (|3.1|) which yields 



2 



1 h 



[ttJ U(l + h/[U/2\) 



(3.2) 
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(where we use |/i(0)| — ► 1 in SC reflecting saturation of the moment). And in 
weak coupling U < 1 (where |/i(0)| = [TJ), §[/%(ft)|/A ~ f/tan-^/t/Ao) is 



likewise insignificant compared to the Zeeman term ft/Ao; although its retention is 
readily shown to be required to recover exactly the leading [/-dependence of the static 
impurity susceptibility Xi(0) (LMA results for which have been considered in ||19| ). 

LMA spectra for arbitrary ft now follow straightforwardly. The zero-field x(0) = 
|[/|//(0)| is first determined using the u = symmetry restoration condition, 
equation (|2.9|) ; and x(h) for any ft > follows as above. £-[■ ft,) (= — cu; ft,)) 
for all is given by equation fl2.8|) with the MF propagators from equation (|2.7|) . 
G a (u;h) = [Q~ x — So-] -1 , G{uj;h) = \'Yl,(jG a (u}]h) and the single-particle spectrum 
.D(u;; ft) then follows directly. 

Before proceeding we comment briefly on the ft-dependence of the spin-flip scale 
u m (h), defined (as for ft = |T6| , |17| ) as the position of the maximum in the transverse 



spin polarization propagator lmll + ~(uj; ft); itself given (§2) by the p-h ladder sum, with 
the bare polarization bubble °n +_ (cG>;ft) expressed in terms of MF propagators. But 
since u m (h) = u m (x(h)) as above, the ft-dependence of the spin-flip scale follows from 
a knowledge of the x-dependence of °n + ~ as considered in [pl|. This may be deduced 
analytically in SC, with the result 

u m (h) =< + 2|/i |ft. (3.3) 

Here |/io| is the MF local moment in zero-field, given explicitly from |/io| = 
2 tan" 1 ^|?7|yU | j for the wide-band AIM [16|; and in practice equation (3.3) is 



numerically accurate for U > 3 or so. In the SC Kondo limit where |/io| 
^m(ft) = + 2ft ~ 2ft for ft/cu^ ^> 1. This is physically correct, it being known 
from solution of the Kondo/s-d model |§ that for fields large compared to the 
Kondo scale u K ~ - and in practice for h/u K > 1 — 10 - the impurity spin becomes 
asymptotically free (albeit with logarithmic corrections to the impurity magnetization 
Mi (ft) H pH|); and for a free spin-| the sole energy scale for spin flips is the Zeeman 
splitting 2ft = — 6j| characteristic of the atomic limit. 

Single-particle spectra for ft = have been considered in references [16-18] to which 
the reader is referred (and are shown as part of figures 2, 5, 7 below). While the 
LMA is perturbatively exact to second order in U in weak coupling |]16| , our primary 
interest is naturally in SC where the Kondo effect prevails. Figure 2a thus shows the 
zero-field spectrum |l!| (solid line) tyA D(uj;0) vs uj/A for U = 6 (wide-band AIM). 
The following spectral features should be noted, (i) The Kondo resonance is correctly 
pinned at the Fermi level, nA D(0; 0) = 1, as follows from self-consistent imposition of 
symmetry restoration (equation Q2.9D ). By the same token its HWHM, the Kondo scale 
cjk, is exponentially small: uk oc uj^ with oc exp(— nU/8A ) in SC as explained in §2. 
(ii) The maximum in the Hubbard satellite(s) occurs at |u;| = j in SC. It is Lorentzian 
in form, but its HWHM is 2Ao and not Ao as simple MF theory would predict, due to 
additional many-body broadening processes (see also below). 
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3.1. Results: h > 0. 

There are three energy scales relevant to the h = ^gfisH dependence of single-particle 
spectra, namely the Kondo scale u^, the hybridization A , and the interaction ~; and 
such that u;^ <C A <C U in SC. There are thus several relevant domains of field strength. 
First the all important Kondo limit, corresponding formally to finite h/u^ and — > 0. 
Here by definition h is vanishingly small compared to A or U; but h/u^ is arbitrary, 
and spans the entire range of field strengths appropriate to the Kondo model. This 
is the universal Kondo scaling regime; we consider it in detail in §4ff. Note however 
that this regime is simply inaccessible to approximate theories in which the Kondo scale 
does not exist (e.g. equation of motion approaches 0); and to those that fail to recover 
an exponentially small Kondo scale in SC, and hence the pristine separation between 
the Kondo scale and Aq or U that is the essence of the Kondo regime (e.g. modified 



perturbation theory fL0|). It is primarily this domain that is of experimental interest in 
the context of quantum dots |3|. || . 

Figure 2 illustrates spectral evolution in the non-universal regimes of field strength. 
For U = 6, and in addition to h = 0, figure 2a shows the LMA 7iA D(u; h) vs u/A for 
V^m = 125 (h/A ~ 0.15 < 1) and h/u^ = 10 3 (h ~ A ); this is continued in figure 2b 
for h/u m = 10 3 , 4 x 10 3 (A < h< §), and h/u^ = 1.2 x 10 4 and 1.6 x 10 4 illustrating 

Consider first the evolution of the 'low'-energy spectral features. By h/u^ = 125 
the zero-field Kondo resonance has split, producing narrow peaks of much reduced 
intensity centred upon |u;| ~ 2h i.e. the Zeeman splitting energy. The latter is also 
found by an equation of motion approach [f|, and seen in DM-NRG calculations by 
Hofstetter |y| for U ~ 3 (see figure 3 of [|lj| for D^(u; h), noting that the field therein 
corresponds to 2h in the present notation and that 2h/A lies in the range 0.1 to 1). 
We add again however that frequencies of order \u\ = 2h lie within the Kondo scaling 
'window' only as 2h/A — > 0, which is not the case for the results shown in figure 2 or 
those of reference |[L5|| . With further increasing field the split peaks associated with the 
erstwhile Kondo resonance remain centred on |a>| ~ 2h, but diminish further in intensity; 
ultimately losing their integrity and being subsumed into the Hubbard satellites as h 
approaches the order U/2 characteristic of the zero-field satellites. This is seen further 
from the LMA ttAqD^(uj\ h) shown in figure 2b inset; it is again qualitatively consistent 



with the DM-NRG results of reference 15 



There are in fact two non-universal regimes of field strength, the crossover between 
which occurs for h ~ ^. This is apparent from the /i-dependence of the Hubbard 
satellites. In SC the latter (figure 2) are centred on |u;| = ¥ + h, as is obvious from the 
atomic limit of the model; and the /i-dependence of which thus becomes significant for 
h ~ C(^) (albeit that satellite shifts are naturally perceptible in figure 2 for h ~ O(A )). 
The position of the Hubbard satellites is however secondary: the /i-dependence of their 
widths is the significant issue, as now explained. 

Upon addition of a |-spin electron to an f-spin occupied impurity, with energy cost 
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Figure 2. LMA spectra nA D(Lj;h) vs w/A for J7 = 6 as discussed in text, (a) 
Upper panel: for h/u>^ — (solid line), 125 (long dash), 10 3 (short dash), (b) Lower 
panel: for h/uj° m = 10 3 (solid), 4 x 10 3 (dotted), 1.2 x 10 4 (short dash), 1.6 x 10 4 (long 
dash. Inset: corresponding irA D^(uj;h) vs w/Aq. 
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~ Y + h corresponding to the position of the upper Hubbard satellite, two subsequent 
hopping processes may occur. The added j-spin may hop off the impurity, itself leading 
to an electron loss rate (and hence HWHM spectral broadening) of Ao- This is elastic 
scattering, and is all that is captured at the one-electron MF level. Alternatively the f- 
spin electron already present may hop off the impurity, leading to a spin-flip with energy 
cost of order u m (h); and hence to an additional loss rate of A provided u m (h) <C j + h. 
This is a many-body process that has no counterpart at MF level. The total loss rate 
of electrons from the site is thus 2A , whence the HWHM of the Hubbard satellites is 
doubled (and their peak intensity correspondingly halved) compared to the simple MF 
result. This behaviour is clearly evident in figure 2 (where ttAqD(u; h) ~ \ at the peak 
positions of the Hubbard satellites); its formal origins within the LMA for h = have 
been discussed in reference [ |16|| . 

For sufficiently large fields however, u m (h) will become comparable to ^ + h and 
the spin-flip energy cost effectively prohibits the additional many-body broadening. 
Since u m (h) ~ 2h for h ^> cu^, this arises for h ~ ^. That the additional many- 
body broadening is 'switched off' for fields of this order is indeed seen in figure 2b for 
h/u^ = 1.2 x 10 4 and 1.6 x 10 4 (h/ [y] ~ 1.5 and 2 respectively); and D(u; h) is given 
asymptotically by 

==E 7 A2 < 3 - 4 > 

2 . [»+'(?+'')] + A » 

such that 7tA /}(cg>; h) | for |a>| = ^ + h. 

Equation ( |3.4p is of course the pure MF spectrum in strong coupling, and that it 
is asymptotically exact for h ^> ^ may be seen from simple consideration of the h- 
dependent one-electron site energies (see equation (|2.1|) ), viz e i(J = ti — ah with = — ~ 
for the symmetric AIM. If ^> and (with the Fermi level) - i.e. if 

/j> \ei\ = y - then in the ground state the impurity is occupied only by |-spin electrons. 



As far as D^(uo; h) is concerned only |-spin electrons are then involved in virtual hopping 
processes; so the fermions are effectively spinless, leaving an effective one-body problem 
with corresponding site energy ef 1 = e, — h. D^(u; h) is then a Lorentzian of width A 



centred onw = £i S (= ~[^+h]), and Di(u; h) follows from p-h symmetry. Equation (|3 
results for -D(o>; h). It is the spectral signature of the free-orbital regime where the static 
impurity susceptibility XiiJ 1 ) coincides with that of the U = limit (as is captured by 
the LMA, see [0). The crossover to such behaviour should be observable in NRG 
calculations of D(uj]h), for although NRG cannot handle adequately the many-body 
broadening of the Hubbard satellites for h ^> ^, one-electron broadening arising from 
the hybridization A(lj) is well captured by a technique introduced recently by Bulla, 
Hewson and Pruschke ||25|| . 

3.2. Approach to the Kondo scaling limit. 

In strong coupling U ^> 1, the low-energy physics of the AIM depends solely upon the 
Kondo scale. The latter itself appears in a variety of superficially different guises, viz 
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Figure 3. ttAqD(u); h) vs uo = Lo/uo m (h) for fixed w m (/i)/w^ = 21 (i.e. h/u^ ~ 10), 
and U = 2 (dotted line), 4 (short dash) and 8 (long dash). The Kondo limit scaling 
spectrum obtained in §5.1 is also shown (solid line). 



))• 

But these are all of course equivalent, being simply proportional to each other (e.g. 



u^, u K or A Z(0) (with Z(0) the zero-field quasiparticle weight, see equation 



A Z(0) = juj^ and luk = 0.691^ within the present LMA |18| ); we denote any of 
them by uo a . 

Although uj a itself depends upon the interaction strength (oj a oc exp(— 7rc7/8A )), 
the fact that it is the sole low-energy scale in SC means that the Kondo/ Abrikosov-Suhl 
resonance exhibits universal scaling in terms of uo/uo a alone, with no explicit dependence 
on the bare material parameters. That the LMA for h = leads to such scaling 
behaviour with progressively increasing U has been shown in reference |16] ; the resultant 



scaling spectrum in the Kondo limit has also been obtained analytically in reference [18 
and shown to give good agreement with h = NRG calculations pOfl . 

Scaling behaviour should likewise arise for h ^ 0, but with the scaling spectrum 
tiA D(lj; h) = F(oj /uj a ]h/uj a ) now dependent upon uj/u a and h/uj a . That such 
behaviour arises within the LMA upon progressively increasing U is illustrated in 
figure 3. We consider a fixed value of u m (h)/u^ = 21, corresponding from equation ( |3.3j ) 
to a fixed h! = h/u^ = 10 in the Kondo limit where | /xo | — > 1 (and to h' within 10% of 
this value for U > 4 - 5). For three different interaction strengths U = 2, 4 and 8, the 
resultant hA D(uj; h) are shown versus u/u m (h) (oc uj/uJ^). They are indeed seen to 
approach asymptotically the Kondo scaling limit (albeit somewhat more slowly than for 
h = [|16|]). The latter is also shown in figure 3, and will now be investigated in detail. 
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4. Scaling spectrum: general considerations. 



Our aim here is to obtain analytically the /i-dependent LMA scaling spectrum 
appropriate to the Kondo limit; and in particular to do so with only minimal assumptions 
about the form of the transverse spin polarization propagator U + ~(uj;h) that enters 
the dynamical part of the self-energy (equation ( p.8|) ). As for the zero-field problem 



considered in reference JL8J the latter will enable us to encompass, and go beyond, the 
particular case considered in the previous section where U + ~(uj] h) is given by the p-h 
ladder sum. 

To obtain the scaling behaviour for the Kondo resonance, one considers finite uj/uj^ 
and finite hi = h/oj^ in the limit oc exp(—nU/8A ) — > 0; the latter projects out the 
non-universal features that are not part of the scaling spectrum (such as the Hubbard 
satellites). Hence, referring to equations (|2.2| , |2.3|) , the 'bare' uj = [lu/u^lu^ = may 
be neglected, as too may the bare field h = h'u® = 0; and A (a;) likewise reduces to 



A(0) = — i sgn(c<j)Ao. The scaling spectrum then follows from equations ( 2.2 



ttA D(lj; h) 



2 ^ 



l + £5?>;/0 



" \h (^(°; h ) + ^ h ) ~ ^ *0l) ] + [ x + A EI >; h) 



as 



(4.1) 



(where t l a (uj- h) = S^w; h) and t^{u; h) - S^(0; h) = S^(cj; h) - S^(0; h) are trivially 
used). Equation ( |4.1|) is general, in the sense that provided the host is metallic it applies 
to any one-electron hybridization A(a>). 

As expected, the scaling behaviour of the spectrum is thus determined exclusively 
by that of the self-energies. And the form of the latter is in turn closely related to the 
corresponding zero-field problem considered in reference ||18|| . Specifically, the transverse 
spin polarization propagator that enters the LMA H a (uj;h) (equation ( |2.8| )) has the 
following functional form in SC, 



with 



-IndT 

7T 



A 



tu m (h) 



mew 



71 



Imn + -(w;/i) = 1 



DC 

Ajdy f(y) 



(4.2) 



(4.3) 



and uj = uj/uj m (h). The essential point here is that Imil + "(w;/i) for h > has the 
same functional form as for h = [18|, scaling in terms of Co = u/uj m (h) in the same 
way as it does in terms of oj/uj^ for h = 0. Such behaviour is physically natural, 
since u m (h) remains the sole low-energy spin-flip scale for /i / 0, just as it is for 
h = where uj^ = uj m (h = 0). Three further points should be noted here. First, 
equation (|4.3|) embodies physically the saturation of the local moment in SC (|//| — ► 1); 
the constant A being h- independent and determined by the resultant 'normalization'. 
Second, by definition of u) m (h), the function f(u) is peaked at uj = 1; and f(u) ~ uj as 
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uj — > 0. Finally, we add that the above behaviour is readily shown to arise explicitly 
with Il + ~(u;; h) given by the p-h ladder sum; from which f(uj) is found to have the form 



m 1* 



UJ 



(4.4) 



1 - 2au + Cj 2 ' 

In the following however, the particular functional form of /(£>) will not be required; we 
shall need it only in §5. 

Using equation ( |4.2| ), the scaling behaviour of E ct (c<j; h) (equation (|2.8|)) follows in 



direct parallel to the h = case [18]. We thus quote only the relevant results from 
reference |T8J. Specifically 



A %(uj; h) = 6(-u)4A I dy f(y) 



(4.5) 



(and £j(co>; h) = E*(— uj; h)), which is required in equation (|4.1| ) and scales solely in terms 
of uj = uj/uj m (h) with no explicit h- or [/-dependence. And for E|(a>; h) (= —Tjf(—u; h)), 



A^Ef^/i) = - In 

7T 



A 






_uj m {h)_ 





A dy f(y)ln\y + uj\ 



(4.6) 



where (see §2) A = min[Z), ¥]. From this the quasiparticle weight, given generally by 
equation (|2~5|) and reducing to Z[K)~ X = — (dT^(uj; h)/du) u= o in the SC scaling regime, 
is thus given by 



1 AA 



and satisfies 



A Z(h) oj m (h) 7r 



Z(h) uj m {h) 



dy 



f(y) 



(4.7) 



Z(0) w° 

(since A is /i-independent). Equation ( ]4.6| ) also yields directly 

oo 

A - x (Sf ( W ; ft) - Ef (0; fe)) = ~A f dy f(y) In 



(4.8) 



1 + - 



(4.9) 



which is likewise required in equation ( |4.1| ) for 7rA Q D(uj; h); and which again scales solely 
in terms of uj. 

The explicit [/-dependence of the zero-field Kondo scale uj^ = uj m (h = 0) follows 
from symmetry restoration equation (^.9|), viz E?(0; 0) = \U (since |/i(0)| — > 1 in SC). 
Equation (|4.6|) for h = then yields directly 



(4.10a) 
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with c aU- and /i-independent constant of order unity given by 



exp 



A / dy f{y)Hv) 



and 



A exp 



-irU 



The remaining quantity required to determine equation ( |4.1| 
A^S^OjA) (equation which reduces in SC to £f(0;/t) 

1 together with 5|/z(/i)| 



£?(0;0), again using |/i(0)| 



Sf(0;/i) 

equation ( ^.2| ) since the bare /i = /i'u;^ = 0). 
then yields A^tf(0; h) = -f ln(a; m (/i)/a;° ); 



(4.106) 



(4.10c) 

for 7tAo-D(c<j; /i) is 
= -f + Sf(0;/i) = 
= (as follows from 



Since A is h- independent, equation 
i. e. using equation 



A^tfi^h) = In 

■ TT 



7T 



Z(0) 



(4.11) 



(with £f(0;/i) = -£f(0;/i)). 

Our final task is to determine the explicit field-dependence of Z(h)/Z(0). This 
may be deduced in two distinct ways. First, and more generally, using an exact Ward 
identity pertaining to the Kondo regime of vanishing charge fluctuations 0: 

dtf (0; h) 



+ 1. 



(4.12) 



dh Z(h) 

Combined with equation (|4.11|) , this leads to a trivial differential equation for Z(h); 
solution of which yields the LMA result for Z(h)/Z(0), viz 

Z(h) , n h 



z(0) 1 + 2M0) (4 ' 13) 

(remembering that we consider finite h/A Z(0) and A ^(0) oc — * 0). Alternatively, 
the /i-dependence may be deduced directly from the LMA with U + ~(ui;h) given 
explicitly by the p-h ladder sum; which we refer to from now on as the LMA(RPA). In 
this case, as discussed in §3, the spin-flip scale u m (h)/u^ a = 1 + 2h/uj^ a in SC. But for 
the LMA(RPA), Af(y) 



5(y — 1) (see QT8J and §5 below), and thus A Z(0) 



7T ,0 



from equation (^4.7|) ; hence using equation (|4.8|) , equation (|4.13|) for Z(h)/Z(0) is again 
recovered. Equation ( [4.13| ) is of course approximate. An exact result for Z(h)/Z(0) 
can however be deduced using the Ward identity equation ( |4.12|) ; and comparison of 
which to equation ( [4.13j ) shows the latter to be a good approximation over essentially 
the entire range of field strengths (figure [TT| ). 

For any finite field h/A Z(0), the scaling behaviour of the single-particle spectrum 
■kAqD{uj] h) (equation 



4])) may now be obtained directly using equations ( |4.5| , 
S7TTL [4. 13| ) together with p-h symmetry. This will be considered explicitly in §5 for a 
particular form of the function f(oj) that determines (see equation ( f4.2| )) the transverse 
spin polarization propagator. First however, we consider predictions arising from the 
preceding analysis that are essentially independent of the details of f{Cj). 
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4-1- Spectral Limits 



We here consider briefly the field-dependence of the single-particle spectrum at the 
Fermi level u = 0, as well as the behaviour of the spectral tails for frequencies 
\u\ = \uj\/uo m (h) > 1. 

The Fermi level spectrum is given from equations (|2.2| , p.3|) by ttAqD(0] h) = 
[1 + ([h — crS^(0; h)]/ Aq) 2 ]^ 1 (with cr£^(0; h) independent of o by p-h symmetry). 
Defining h = h/AoZ(0), the LMA result for 7rAo-D(0; h) in the Kondo limit of finite h 
and A o Z(0) — > thus follow from equations Q4.11J , |4.13| ) as 



%A D(0;h) 



1 



4 



In 



7T 



2 . 



-i 



with asymptotic field dependencies: 

7iA D(0;h) I -Ah 2 



7iA Q D(0;h) 



hy$>i 



'4 

- ln(/i 

7T 



-2 



(4.14) 



(4.15a) 



(4.156) 



The exact behaviour of 7rA -D(0; h) can be obtained, since the (excess) impurity 
magnetization Mi(h) follows generally from the Friedel sum rule as n M\(h) / g /ib = 
tan- 1 [(/i - ctS^(0; h))/A }; whence 



nA D(0]h) = cos 2 



7rMj(fe) 

g^B 



(4.16) 



with Mi(h) known from the Bethe Ansatz (BA) solution of the Kondo/s-d model EI 



From this the exact low- and high-field asymptotics of ttAqD(0; h) are readily obtained, 
and are given by equations (4.15). The LMA thus captures these correctly, which is 
non-trivial. 

A full comparison between the LMA and exact BA results for 7rA -D(0; h) is given 
in figure |], and the level of agreement for all field strengths is self-evident. We also 
add that numerical renormalization group calculations of the single-particle spectrum 



for the Kondo model [14] yield excellent agreement with the BA result for irA D(0; h), 
over the modest range of field strengths considered in reference [Q. Static properties 
arising from the LMA are discussed further in reference | I9| . In particular the Wilson 
ratio Rw(h) — 2 V h M, |22fl is correctly recovered; and the field dependences of M^h) 



and the corresponding static susceptibility Xi(/i) are likewise found to be asymptotically 
exact in both the weak and strong field limits. 

The behaviour of the spectrum for frequencies \u\ = \u\/u m (h) ^> 1 is also 
readily obtained. From equations (|4.5| , f4.2| ), Aq 1 Sj(cj; /i) is given asymptotically by 
A^T^iuj; h) = 40(-o>). Equation Q for \u>\ > 1 yields A : (Sf (u; h) - (0; h)) = 
~ ^ Io° f (u) ^ n [\ti\ / y] = — -ln[|a)|c] where equation ( |4.3| ) is used and the constant c 
is given by equation ( [4.10£|) ; and Aq 1 E| l (0; h) = — ~ \n[uj m (h) /cu^} from equations ( |4.11| , 
Hence Ef{u;h) = Ef{0;h) + [Ef{u;h) - Ef (0; h)} is given for \uj\ > 1 by 
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Ao 1 



Ef(u,;/0 = -iln(| 



where equation ( |4.10a| ) is used, and the [/-dependence 



of cu' m is given explicitly by equation ( |4.10c| ) 
follows from equation (|4.1| 



The scaling spectrum for \u\ ^> 1 then 



as 



7iA D(u] h) 



H>1 



+ 



+ 25 



(4.17) 



[J In Mr + 1 [J IBM]' 

where a/ = uj/u' m ; equivalently, the frequency dependence may be recast in terms of 
uj/uj k (with uo K the HWHM of D(uj; 0)), using u K /u' m = 0.691 found for the h = LMA 

El- 

Equation ( [4.17 ) is precisely the result for the behaviour of the scaling spectrum 
'tails' obtained in reference |1§ for h — 0, and there formally applicable for |u;| 3> 
co>^ = u) m (h = 0). The resultant slowly varying tail is logarithmic in form, as opposed 
to the algebraic Doniach-Sunjic |26| decay D(u;0) ~ (\uj\/uk)~2 hitherto thought to 
arise 



(see e.g. [g7| Q). We have however argued |L8| that the logarithmic decay is 
entirely natural in physical terms, and shown both that equation ( |4.17|) gives excellent 
agreement with NRG calculations of D(u;0) (see e.g. figure 2 of reference |18[); and 
that this long tail in fact dominates the h = scaling spectrum, the crossover to Fermi 
liquid form occurring only on the lowest energy scales |o;|/u;k *C 1. The arguments 
leading to equation ( 4.17] ) show additionally that this tail behaviour is asymptotically 
common both to h = and h > 0, where it arises for |d)| = \u}\/u xa {h) ^> 1; i.e. for 
^> 1 + \h (using equations ( |4.S| , f4.11| )). This will be seen explicitly in the results 
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shown in the following section. 



5. LMA Scaling spectrum. 



To obtain the LMA scaling spectrum on all energy scales requires a full specification of 
the transverse spin polarization propagator (equation [4.2|) . This is now considered, first 
(§5.1) using the LMA(RPA) with U + ~(lj; h) given explicitly by the p-h ladder sum; and 



then (§5.2) via a simple modification thereof introduced previously in reference [jig] . In 



}5.3 we consider the spinon approximation to the /i-dependent single-particle spectra 



developed in reference JTT|. Comparison to recent transport experiments on quantum 



dots in the Kondo regime || is made in §5.4. 



5.1. LMA(RPA). 



For the LMA(RPA) the function f{Cj) that determines ImIT + (u;; h) (equation ( [4.2|) ) 
has the form equation ( |4.4|) , where the x-dependence of a (and A, see equation ( |4 . 3|) ) 



is given explicitly in reference [jlqj . From this it is known that in SC U 3> 1, 
a — > 1 and A ~ [2(1 — a)] i/ir — > such that Af(y) = 5{y — 1) i.e. from equation ( |4.2| ) 
ilmn + ~(c<j; h) = 5{uj — u m (h)) reduces to a delta function centred on u m (h) = uj^ + 2h. 
Note also that u® a = u' m follows from equations ( f4.1Qq| , f4.10£| ), with the [/-dependence 
of u' m given explicitly by equation ([4.10c]) . 



From equations ( |4.5| , |4.9|) it follows directly that 

Aq Hzfiu; h) - S?(0; h)) = -- In \u + 1| (5.1a) 

7T 

Ao^K^j/i) =40(-[w + l]) (5.16) 

where uj = u/u m (h); while Aq 1 E| 1 (0; h) = — ^ \n[u m (h) / u>^} from equations ( f4.11| , [4 . 8|) . 
The spectrum in closed form then follows simply from equation ( [4.1|) ; specifically for 
u; > (since D(—u; h) = D{uo\ h)) by: 



itA D(uj; h) = 

if 1 + 1 + 46(u> - u> m (h)) 

2 \ [i In |a/ + <(/0 1] 2 + 1 [J In |a/ - <(/,) |] 2 + [1 + 40(jS - u' m (h))] 2 

Here to' = uj/ui' m , and likewise uj' m (h) = u m (h)/uj' m = 1 + lb! where b! = h/uj' m (= h/u^ 
as used in §3); equivalently oj' m (K) = 1 + jh, where (as in §4) h = h/A o Z(0) and 
A Z(0) = juj' m for the LMA(RPA). Equation ( p.2| ) is the Kondo limit scaling spectrum 
shown in figure ^| above (for fixed w' m (h) = 21) and arising as the large-?/ limit of the 
AIM. It naturally reduces in the zero-field limit (where w' m (0) = 1) to the result obtained 
previously in reference |TS[ . 

The field dependence of the resultant spectrum is shown in figure || nA D(u; h) 
vs u/luk for h = h/A o Z(0) = 0, 0.2, 0.5, 1, 2.5 and 5; here as throughout, wk is 
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Figure 5. LMA(RPA) scaling spectrum itA D(uj; h) vs uj/ojk for h — h/A n Z(0) = 0, 
0.2, 0.5, 1, 2.5 and 5 (from top to bottom). 



defined as the HWHM of the zero-field spectrum (with ujj^/u' m = 0.691 for the LMA 
pl|). The spectrum at the Fermi level, u = 0, decreases monotonically with increasing 
field (figure |J); and for sufficiently small fields {d 2 D{uj; h) / du 2 )^^ < 0, whence the 
maximum in the Kondo resonance remains at u — 0. As pointed out by Costi in a 
recent NRG study of the Kondo model irjj however, the Kondo resonance 'splits' above 
a certain field H c where the us = curvature changes sign. For the LMA(RPA) this may 
be obtained analytically, and occurs at 2h c /uK = g^BHc/u^ = 0.600; compared to the 
corresponding NRG value [14| of 0.5. 

With increasing h > h c the split Kondo peaks in D(uj; h) move further apart, 
although as seen (figure H) the generic tail behaviour equation (|4.17|) is always 
approached at sufficiently high frequencies. The LMA(RPA) peak splitting is readily 
deduced from the individual D a (u;h) (with tiA Di(uj; h) given in its entirety by the 
second term in equation (|5.2|)): since uj' m {h) = 1 + 2h', it follows directly using 
equation (|5.2|) that D a (u]h) = D a {oj + 2ah; 0) with D a (u;0) the zero-field case. 
D a (u; h) thus amounts simply to a rigid shift of its zero-field counterpart, and hence 
has its peak maximum at u' p = —2ah' precisely. The corresponding maxima in D(u; h) 
(= \ J2 a D a (oj] h)) for h' > h' c likewise occur in practice very close to \u p \ = 2h', figure ^| 
(and asymptotically so for h! 3> 1); with a total peak splitting of Ah, or twice the Zeeman 
energy. We add that a peak position of \u p \ = 2h in D a (u; h) arises as the large-field 



(h ^> luk) asymptotic behaviour within the spinon approximation [p t| (discussed further 
in §5.3). It also arises within an equation of motion approach 0; although the apparent 
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parallel here is somewhat misleading since the Kondo scale does not exist within such 
an approach, and hence neither does the Kondo scaling regime of finite h/u>K and uj/uj-^. 

5.2. Beyond the LMA(RPA). 

While the LMA(RPA) captures well the spectrum at the Fermi level (figure f|), and 
probably also the spectral tails (as for h = [)18| , see also figure ^| below) , its quantitative 



limitations are apparent from figure ||. The split Kondo peaks, for example, are too 
pronounced in comparison to recent NRG calculations for the Kondo model fL4| . The 
origins of this limitation stem from the divergence in the LMA(RPA) h) at \oj\ = 1 

(equation ( 5.1cj) ), reflected by the 'dip' in 7tA D(uj; h) at uj — u> m (h) evident in figure |. 



As known for the h = problem [16|, [Tj|, this is entirely an artifact of the specific 
RPA-like form for the polarization propagator, embodied in the fact that the resultant 
Imil +_ (ti;; h) is a 5-function at uj = uj m (h): in reality one expects ImlT (a;; h) to have 
a finite width. 

To rectify this deficiency we proceed as in reference []18|| . We retain the form 
equation (f4.4|) for /(d)), which has a finite width provided q^I; and we employ a high- 



frequency cutoff uj c to render f(oj) normalizable (equation ( [4. 3D ). The width parameter 
a is then determined by requiring that the leading low-frequency behaviour (oc uj 2 ) of 
the imaginary part of the conventional single self-energy for h = is recovered exactly. 
As discussed in reference JL8J this requires A = |[^/ A Z (0)] 2 (= |[o; m (/i) / A Z (h)} 2 , 



see equation (fig) ), which via equations ( f473| , |4~7| ) in turn implies a simple equation that 
determines a uniquely for the chosen cutoff u c . The latter is of course arbitary but, as 



expected physically, results are not sensitive to it |L8|]; in practice, as in reference [18 
we choose uj c = 10 (and hence a = 0.308). 

While the effects of this modification are rather minor in comparison to the 
LMA(RPA) forfe = pfl, they are more significant at finite fields. For h = h/A Z(0) = 
1, figure |6| shows hA D(oj; h) vs uj/ujk for the LMA(RPA) compared to the LMA outlined 
above (and refered to simply as the LMA). The zero- frequency and asymptotic tail 
behaviours of the two naturally coincide, being independent of the detailed form of 
f(oj) (see §4.1). But elsewhere the differences are clearly quite significant. In particular 
the Kondo peaks in the LMA are less pronounced, and for the example shown in figure ^| 
their splitting is smaller in comparison to the LMA(RPA). 

In figure the /i-dependence of the resultant LMA spectrum is shown (for the same 
fields as in figure | for the LMA(RPA)). For h < 1, the first effect of the field is to erode 
the zero-field Kondo resonance 'on the spot' - first diminishing and then splitting the 
resonance, but doing so largely under the envelope of the zero-field spectrum itself. With 
further increasing field however, and again in contrast to the LMA(RPA) (figure |), the 
split Kondo peaks broaden and diminish further in intensity; and their maxima move 
outside the zero-field spectral envelope, although the tail behaviour of the resonance is 
asymptotically common for all fields (as shown generally in §4.1). 

The LMA spectra shown in figure [7] agree rather well with those obtained 
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Figure 6. For h = 1, comparison between scaling spectra irA D(uj- : h) arising from 
the LMA(RPA) (dashed line) and the LMA (solid line) discussed in text. 
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Figure 7. LMA irA D(uj;h) vs lu/uj k for h = h/A o Z(0) = 0, 0.2, 0.5, 1, 2.5 and 5 
(from top to bottom). 
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from recent NRG calculations |TJ[ (see figure 3 therein). In particular the spectral 
characteristics outlined above are as found in the NRG calculations, and the qualitative 
similarity of figure [7] and figure 3 of reference [13] is self-evident. Detailed comparison 



may be made by noting that the LMA spectra at uo = for a given value of h = h/A Z(0) 
(given explicitly by equation ( f4.14| )), are in very good quantitative agreement with the 



uj = NRG spectra at the same value of the ratio H/T^ employed in reference \T%\. The 
resultant comparison as a function of frequency is not of course quantitatively perfect; 
but it is certainly rather good, and the LMA is to our knowledge the only theoretical 
approach that bears such comparison to the NRG data (see also figure || below). 

With increasing field the split Kondo resonance becomes increasingly broad/diffuse 
(see figure [7| and figure pi] below), and we now comment on the large- h dependence of 
the peak maxima in the LMA D(u; h) (or equivalently in the -D CT (u;; h)). This we find 
numerically to be of the form \lv p \ ~ hlnh for h = h/AoZ(0) ^> 1. Hence |o; p | exceeds 
the value of 2h (the Zeeman splitting energy), generally regarded as the maximum 
spectral shift: see e.g. HIT , 14| . For the physical reasons outlined in §3, it is we believe 



clear that the maximum shift in the transverse spin polarization propagator - which 
directly probes spin-flip dynamics - will not exceed an energy on the order of the free 
spin-i Zeeman splitting 2h. We do not however know of a convincing reason, physical or 
otherwise, as to why this should also hold for the single-particle spectrum in the Kondo 
scaling regime (note that the NRG calculations of reference [H|, which report |u; p | < 2h, 



are limited to comparatively small fields of H/T^ < 5 or so). This issue should however 
be resolvable by NRG calculations, possibly necessitating use of the DM-NRG technique 
introduced recently by Hofstetter flT5| |. Indeed we add that preliminary calculations on 



the AIM in the strong coupling scaling regime, using the DM-NRG method, confirm 
spectral shifts in excess of 2h for sufficiently large fields [^3|, as seen directly in 
figures H, HO below. 



In the opposite limit of h — > 0, an exact result for the spectral maximum 
u> p (h) = — a\u> p (h) \ of the cr-spin spectrum D a {uj; h) may be obtained using Fermi liquid 
theory, as discussed in the Appendix. This leads to 

M = ITpEZW h (5 ' 3) 

where (3 = lim/ w0 [hm a; _ ) o(Iv(u ; ; K)/uj 2 )]. In the trivial non-interacting limit where 
-Rw(O) = 1 (and (3 = 0), this recovers correctly |u;p(/i)| — h (— ^gfiBH). In the 
strong coupling Kondo limit by contrast, i?w(0) = 2 and the exact (3 = l/(2Ao^(0) 2 ) 
0; hence 

\u p (h)\ = | h. (5.4) 



Equation (|5.4|), valid for h <C 1, is exact for the Kondo model (albeit somewhat 



unexpected, conventional lore suggesting |cG> p (/i)| = 2h). It is not however recovered by 
the LMA for D a (uj; h) which instead yields |d;p(/i)| = h as h — > 0; reflecting the fact that 



while i?w(0) = 2 arises correctly within the LMA JT9|, the resultant (3 is readily shown 
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Figure 8. LMA 7rA Dj.(a;; h) vs cj/^k for h 
(point dash); and itA D(lu;0) (solid line). 



1 (short dash), 2.5 (long dash), 5 



to be twice the exact value. Neither does equation ( f5.4| ) appear to be recovered by 
the spinon approximation to the single-particle spectrum [|TT] (discussed further in §5.3 
below), the field dependence of \u p (h)\/2h being shown in figure 2 (top) of reference (Tl[ . 
For the lowest field considered there {q^H/Tq = 2h/uj^ ~ 1), \u p (h)\/2h is close to 
but slightly less than |, and appears to be diminishing further with decreasing field 
(although its — ^ asymptote has not to our knowledge been determined). NRG 
calculations of the Kondo model were originally reported |14] to yield |a> p (/i)| ~ 2h for 
h <C 1, but we understand that subsequent reanalysis of the low-field data J30[ is now 
consistent with the exact result equation (|5.4T) . 

Finally, recall from §2 that for h = 0, D(u;0) given from equation (|2.13|) is 
independent of spin a, and is equivalently the zero-field a-spin spectrum. That the 
present LMA does not recover correctly the above low-field spectral shift of D a (Lv; h), 
stems from the fact that the individual LMA D a (uj;h) as h — > do not separately 
coincide with D(u; 0) for arbitary u (although the differences in general are rather minor, 
and -Do-(0; h) = D(0; h) for all h and either a). This in turn is a natural consequence 
of the fact that (§2), for any h ^ 0, one or other MF state is picked out according 
to sgn(/i); and as a result the peak maxima in the individual D a {uj; h) for small h are 
not captured correctly. Note that the above remarks apply solely to the individual 
D a (u; h), since D(u; h) for all u evolves continuously in h to its h — limit D(uj] 0) (see 
§2). It is in fact possible to recover the correct low- field shift in D a (oj\ h) within an LMA 
framework; but discussion of this would take us too far afield, and the present LMA 
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Figure 9. 7tAqDi(uj; h) vs lo/u>k for /i/wk = 25. Solid line: LMA. Dashed line: result 
from DM-NRG calculations lEl. 



does appear otherwise to account rather well for the overall spectral behaviour of the 
individual D a (u; h). These are illustrated in figure ^| for h = 1, 2.5 and 5, and exhibit 



the same behaviour as the NRG results shown in figure 2 of reference |14| . In addition, 
the LMA for ttAqD^co; h) is compared explicitly in figure |] with results arising from 
the DM-NRG approach [|15| for a field Ji/uk = 25 (kindly provided by W. Hofstetter 
23| , and obtained for the AIM at a strong coupling U = 8). The LMA is seen to 



account well for the DM-NRG data. Note in particular that both approaches clearly 
yield peak maxima \u p \ that are in excess of the Zeeman splitting 2h; and that, on the 
low-frequency side in particular, the slow logarithmic tails in 7rA .D|(a;; h) (given by the 
first term on the right hand side of equation ( 4.171) ) are clearly evident. 



5.3. Spinon approximation. 

Here we compare the LMA results of the preceding section to the spinon approximation 
(SA) JO]], in which the /i-dependent single-particle spectrum is approximated by the 
density of states for spinon excitations obtained via the Bethe ansatz; and make 



further observations on the SA itself, based upon the results presented in reference [II 



We denote the SA to the single-particle spectrum by p a iy>', h) (in the notation of 
reference fill ); with p a (u; h) oc D a {uj; h). 

To compare directly the LMA and SA we require the full proportionality, p a (u); h) = 
CD a (u] h) with C an uj- and /i-independent constant. This is clearly given by 
C = p (T (0;0)/D ff (0;0) = p (T (0;0)7rA (using ttAqA^O; 0) = 1). It is however known 
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1 3lfl that the zero- field spinon spectrum p a (u; 0) is a pure Lorentzian, with a HWHM 
denoted by T in reference |TI[ (such that T = wk in the notation of the present work). 
Hence p a (0;0) = l/vrT , and C = A /T . p(u>;h) = ^2 a Pa(u;h) is thus related to 
D{u;h) {=\T^D a {u- 1 h))hy 

p(uj;h)2h = -— nA D(uj;h). (5.5) 

7T Co"k 



In figure 1 (bottom) of reference fll |, and for a field 2h/T = 48 (= 2/i/o;k), 



p(uj;h)2h vs u/2h is shown in the interval \uj\/2h < 1 (to which frequency range the 
SA is restricted in practice |TT|). In figure [K] this SA result is compared to that arising 
from the LMA (via equation ( |5.5| )). The two clearly differ very significantly for all u. 

To gain some insight into the above, it is natural to consider the functional form of 
the SA Pj(a;; h) (such that p(u; h) = Pi(ui; h) + (— u>; h) by p-h symmetry). In figure 1 
(top) of reference |TT|, pi(u>] h)2h vs u/2h is shown, again for h/T = 24. Although not 
remarked upon in reference ]TTJ, p^(oJ\ h) appears to be a Lorentzian, given by 



(5.6) 



where (using the notation of reference |IT[| ) E mSLX (h) and AE(h) denote respectively the 
peak position and FWHM. The numerical validity of equation (|5.6j ) may be confirmed 
by direct transcription of the data given in figure 1 of reference |l ij for h/T = 24; the 
Lorentzian form equation ( |5.6|) fits the data highly accurately (the correlation coefficient 
being 0.99995). 

That p a (u>; h) should be a Lorentzian is not perhaps surprising, it being well known 
]3jfl that the zero-field p a (u;0) is a pure Lorentzian; and although verified directly 
only for h/T = 24, we naturally assume p a (u>; h) to be a Lorenztian for general h/T Q . 
Granted this, the SA to the uj = single-particle spectrum -D(0; h) (= D a (0] h)) then 
follows directly from equations ( |5.5| , |5.6| ) as: 

T (AE{h)/4h) 



irA D(0;h) 



2h 



t (h) 



2i, 



+ 



AE(h) 
4h 



(SA) 



(5.7) 



The /i-dependence of D(0;h) is however known exactly (see equations ( [4.16| , 4.15)), 
enabling the validity of the SA to be ascertained from the /i-dependence of E m3X {h) and 
AE(h) given in reference | i~If . Defining h = h/T , their asymptotic behaviours are 
known analytically for h ^> 1 in particular [[□]]: E may[ (h)/2h ~ 1 — (9([ln(/i )] _1 ) and 
AE(h)/2h ~ (^([ln^/io)]- 1 ). Hence from equation the SA to tiA D{0; h) for large 
fields ho ^> 1 is 7rAo-D(0; h) ~ [/io ln 2 (/io)]~ 1 . This is not however the exact asymptotic 
behaviour, which from equation ( [4.15£|) is in contrast given by ttAqD(0; h) ~ [ln 2 (/io)] 1 
(using ho oc h = h / A Z (0)) . And this in turn underlies the disparity evident in figure |10| 
between the SA and LMA results for D(0; h) in particular (for fields of this order the 
LMA is very close to the exact -D(0; h), see e.g. figure £|). 
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w/2h 



Figure 10. Spinon approximation to the single-particle spectrum jTl]], p(u>\ h)2h vs 
uj/2h (dotted line), compared to the corresponding LMA result [Ah/ ■kuj-k\itIS.qD{uj] h) 
(solid line); for a field h/u>K = 24. The experimental data of reference || is also shown 
(dashed line); see text (§5.3) for full discussion. Inset: LMA spectrum on an expanded 
frequency scale (solid line), compared to [Ah/ ttu)yl\it&.qD(lu\ h) vs u>/2h obtained via 
the DM-NRG method H |§ for h/u) K = 25 (dashed line). 



The deficiencies of the SA are apparent even for h = 0, where the resultant 
Lorentzian spectrum WjR is simplistic, failing to recover the logarithmic tails known 



to dominate the single-particle scaling spectrum (see e.g. [18]). The approximation 
does not moreover appear to improve for h ^ 0, as evident from the discussion above. 
In this regard it is perhaps salient to note that the large-field asymptotic behaviour 
~ [h ln 2 (/io)] -1 deduced above for the to = spinon spectrum, is in fact that of the 
exact excess impurity susceptibility XiW 01 • But whatever its physical content, the 
SA does not in our view provide a qualitatively satisfactory description of the single 
particle spectrum. 

5.4- Experiment. 

The differential conductance of a quantum dot in the presence of a magnetic 
field, G C (V; h) with V the applied (drain-source) voltage, has recently been studied 
experimentally by Goldhaber-Gordon et al 0]. Their data for an applied field H = 7.5T 
has been compared to the SA results JTTJ, taking gfi B H/u>K = 2h/uK = 48 and \g\ = 0.36 



(as suggested by ESR measurements on 2DEGs [^j). With this we likewise compare 
the experimental data |§ to the LMA results of §5.2. The two steps involved are as 
follows. First, considering T = 0, we take (e 2 /27rh)~ 1 G c (V] h) oc nA D(uj = V;h). 
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This proportionality holds strictly for the linear differential conductance (V = 0), non- 
equilibium effects arising for V 7^ thereby being neglected in practice (as also in e.g. 



reference [if]). Second, the proportionality constant 7 is determined simply fom the 
experimental linear differential conductance, viz 7 = [(e 2 /2nh)~ 1 G c (0; h)]/(7rA D(0; h)) 
with (e 2 /27rh)~ 1 G c (0; h) taken directly from the data in reference 0] and D(0;h) the 
LMA spectrum at u = 0. 

The latter step requires an explanation, since the LMA spectrum is for T = 
while the experiment |§ is performed at T = 90mK (approximately twice the Kondo 
temperature). The Kondo resonance is continuously destroyed by the separate effects of 
temperature (as controlled by the ratio T/uk) and an applied magnetic field (controlled 
by gfi^H/uj^). But if g^H ^> T the latter effect overwhelms the former, and 
comparison to the T = limit is thus warranted (the validity of which argument is 
in fact consistent with the NRG results of reference |I4[], figure 3). This is indeed the 
relevant case in the experiment at H = 7.5T and T = 90mK, where (with \g\ = 0.36 
as above) g^H/T ~ 20. 



The resultant comparison between the LMA and experiment H] is shown in figure f 



as a function of u/2h (with [4h/7ru)K\irAoD(u); h) shown to enable comparison to the 



SA results [^TJ also given in the figure). The agreement with experiment is on the 
whole rather good, particularly for |u;| < h where the uj (= V) dependence of the 
experimental G C (V; h) appears to be quite well accounted for by the LMA. Insofar as 
the modest differences between experiment and the LMA may be attributed to non- 
equilibrium effects, neither these nor (as argued above) thermal effects would thus 
appear to dominate the experimental observations. In this regard our interpretation 



of experiment differs significantly from that of reference [If ; which is natural given the 



marked disparity (figure [Tt|) between experiment and the SA results [|TTJ. The LMA is 
in turn in rather good agreement with results obtained from the DM-NRG approach 



Tq , |23| , particularly for |u;| < h. This is seen from the inset to figure [10] where the 



above LMA results are compared to DM-NRG data for 2h/uj^ = 50 (see also figure §). 
6. Summary. 

The subject of this paper has been single-particle dynamics of the symmetric Anderson 
impurity model in the presence of a magnetic field. While topical, the problem certainly 
poses a significant theoretical challenge. The non-crossing approximation [33-35,24] for 
example, which despite its inability to recover Fermi liquid behaviour at low-energies 
has been used to considerable effect for H = (see e.g. 0), fails quite dramatically for 
H 7^ 0; producing spurious spectral peaks at the Fermi level that are symptomatic 
fTf| of its origins as a large- N theory. Likewise equation of motion approaches simply 
lack the low-energy Kondo scale, modified perturbation theory |TD[ is confined to weak 



coupling, and the spinon approximation [ 1 I ] suffers from the qualitative limitations 
discussed in §5.3. 

The local moment approach [16-19] developed here transcends many of the 
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limitations of previous theoretical approaches. All energy scales, field and interaction 
strengths are handled by it, leading thereby to a rather comprehensive description of 
the problem that recovers in particular the important strong coupling, Kondo scaling 
regime. Here, as for the zero-field case [JIB] , the LMA appears to pass the acid test 



of comparison to benchmark numerical calculations provided by the NRG JTJ], [T5|, [23|] ; 
as well as yielding rather good agreement with experiments on quantum dots [0]. Its 
strengths stem in part from its inherent simplicity and physical transparency, together 
with the fact that it is not confined e.g. to problems that are ubiquitously Fermi 
liquid-like on low-energy scales |L7|, ^Ofl . As such, we believe it provides a powerful tool 



for further investigation of a wide spectrum of quantum impurity physics, and related 



lattice-based models within the framework of dynamical mean- field theory |36], p7 | 
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Appendix 

Using microscopic Fermi liquid theory, we obtain two exact results for the Kondo limit 
of the AIM; specifically for the field dependence of the quasiparticle weight (Z(h)/Z(0)) 
and the asymptotic low-field behaviour of the spectral shifts (uj p ) in D a {uj\ h). 

The basic underlying equations are as follows. The (excess) impurity magnetization 
Mi(h) follows directly from the Friedel sum rule (see e.g. @); and for the symmetric 
AIM is given generally by 

Mm = — tan" 1 \{h - o-l£(0; h))/A ] (A.l) 

(where <tS^(0; h) is independent of spin, see equation fl2.4|) ). From this the 
corresponding impurity susceptibility follows, Xi(^) = \9t>>~B{dMi/ dh). The Wilson 
ratio Rw(h) = cxi(h) / '^(h) , with 7i(/i) the linear specific heat coefficient (c = 
[27Tk B } 2 /[3(gfi B ) 2 }); it is given by 

Rw(h) = Z(h)[l - o-(dtf(0- h)/dh)}. (A.2) 

And using ( |A.1| , |A.2|) , Xi(h) may be expressed as 

Xi{h) = — - — £(0; h) (A.3) 

with D(0; h) = D a (0; h) the spectrum at the Fermi level u = 0. Finally, in the Kondo 
regime of vanishing charge susceptibility, we will employ the Ward identity || 



1 | T f d^(0;h) \ _ 2 f dE*(u;h) 
dh I \ du 



(A.4) 

LU=0 
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Figure 11. Exact field dependence of the quasiparticle weight, [Z(h)/Z(0)]/(l + ^h) 
vs h = h/AoZ(0) (solid line); the large-/i asymptote (equation ( A. 9 )) is also shown 
(dashed line). 



(which is not specific to the symmetric model). 

Consider first the Wilson ratio. Since Z(h) is given by equation ( |2.5| ), 
equations ( |A.2| , |A.4| ) yield directly Rw(h) — 2 V h. For h = 0, this result is well 
known since the classic work of Nozieres ]?]]. Its validity for all h was first suggested by 
Wiegmann and Finkelstein p2 1 and later proven using the Bethe ansatz || . It is not of 
course specific to the symmetric model: the generalization of (|A.2|) for the asymmetric 
AIM is readily shown to be R w (h) = £ ff D a (0; h) [l-<r(0f£(O; h)/dh)]/ £ CT .0,(0; h)[l- 
(dT^(u; h) / du) u= o\; combined with (|A.4j) , Rw(h) = 2 V h again results. 

2, 



We now consider Z(h)/Z(0). From (|A.2|) , using Rw(h) 
Z(h) [l-a(d£*(0;h)/dh) h = o ] 



Z(0) 



[l-a(d£*(0;h)/dh)} 



(A.5) 



With £*(0; h) given from flftlp , 

1 - a{dtf{0; h)/dh) = sec 2 

and hence 



nMi(h) 



2ttAo , 



(A.6) 



Z(h) 
Z(0) 



cos 



7rAfj(fe) 
9MB 



(A.7) 
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Here Xi(h) = Xi(h) /Xi(fy> an d from 
Xi(O) = ( gf i B ) 2 /{irA Z(0)) (= (gfi B ) 2 /4kT L where kT L 



using i? w (0) = 2 and 7rA £>(0;0) = 1, 
f A o Z(0)). Equation (TO]) is 



exact, and M\(h) (and hence ls likewise known exactly from the BA solution of 

the Kondo model PH| . The h = h/A Z(0) dependence of Z(h)/Z(0) may thus be found 
explicitly, and is shown in figure [11] where [Z(h)/Z(0)}/[1 + %h] is plotted; this being 
the ratio of the exact Z(h)/Z(0) to that arising within the LMA (see equation ( [4.13] )). 
The latter is thereby seen to concur well with the exact result, deviation from which is 
typically < 10% save for h ~ 0(1) where it is slightly more. 

The large-/i behaviour of Z(h)/Z(0) bears note, it being known from the BA 



solution that [21 



Mi(h) 



1 - 



1 lnln(foA) 
2\n(bh) 4[ln{bh)] 2 



+ 0([\nh] 



(A. 



(where b = A^fe/n). The LMA for M^h) (discussed in |TJ) recovers correctly 
the leading logarithmic approach to saturation ([2 ln(/i)]~ 1 ), but not the subleading 
In ln(/i)/[ln(/i)] 2 corrections. It is in fact the latter that generate the leading logarithmic 
corrections to Z(h)/Z(0) ~ for h ^> 1; which using (|A.8|) and ( |A.7| ) are given by 

Z(h) 



/»>i 



Z(0) 



2 



1 + 



21n(6A) 



(A.9) 



The asymptotic behaviour ([A.9Q is also shown in figure |TT| (dashed line). 

Finally, we obtain an exact result for the asymptotic low-field behaviour of the 
spectral shifts. G a (uj] h) itself is given by equation ( |2.3|) , with S tT (a;; h) the exact self- 
energy; and we consider explicitly the wide-band AIM for which A(u) = — i sgn(o;)Ao 
(this does not of course impose any restrictions on results for the Kondo regime). 
D a {uj]h) = — 7r _1 sgn(u;)ImG' -(u;; h) follows from equation ( |2.3|) , and the spectral 
maximum therein arises for uj = u p (h) such that (dD a (uj; h)/duj) u = 0; to which 
equation we seek the asymptotic low-field solution with u p oc —ah, say uj p = —jah. 
This is merely a matter of algebra: one needs simply a standard low-frequency expansion 

of the self-energy, £*(w; h) - S^(0; h) [Z(h)- 1 - l]u and t\(u; h) ~ (3(h)uj 2 ; and 

to recognize from QS^) that lim^ [l - f S^(0; h))Z(h) = R w (0). With (3 = (3(h = 0), 
the general result for the AIM is thereby found to be 

, , Rw(0) 



[l + (3A Z(0y 



rh 



(A.10) 



as discussed in 55.2. 
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